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To help understand the high activity of silver as an oxidation catalyst, e.g., for the oxidation 
of ethylene to epoxide and the dehydrogenation of methanol to formaldehyde, the interaction and 
stability of many different oxygen species at the Ag(lll) surface has been studied for a wide range 
of coverages. Through calculation of the free energy, as obtained from density-functional theory and 
taking into account the temperature and pressure via the oxygen chemical potential, we obtain the 
phase diagram of 0/Ag(lll). Our results reveal that a thin surface-oxide structure is most stable 
for the temperature and pressure range of ethylene epoxidation and we propose it (and possibly 
other similar structures) contains the species actuating the catalysis. For higher temperatures, low 
coverages of chemisorbed oxygen are most stable, which could also play a role in oxidation reactions. 
For temperatures greater than about 775 K there are no stable oxygen species, except for the 
possibility of O atoms adsorbed at under-coordinated surface sites (i.e., imperfections, defects). At 
low temperatures (< 400 K at atmospheric pressure), provided kinetic limitations can be overcome, 
thicker oxide-like structures are predicted. Due to their low thermal stability, however, they can be 
ruled out as playing an important role in the heterogeneous reactions under technical conditions. 
Bulk dissolved oxygen and a molecular ozone-like species adsorbed at a surface vacancy, as have 
been proposed in the literature, are found to be energetically unfavorable. The employed theoretical 
approach for calculating free energies and predicting the lowest energy structures in contact with 
species in the environment ("afe initio, atomistic thermodynamics"), affords investigation of a system 
that seemlessly connects standard (T — K) DFT results, characteristic of "typical" theoretical 
surface science studies, through to those valid for the conditions of catalysis. Though the error bar 
of the noted theoretical temperatures is noticeable (± ~55 K), the identified trends and physical 
descriptions are useful. 

PACS numbers: PACS: 82.65.Mq, 68.35.Md, 68.43.Bc 



I. INTRODUCTION 

The identification of active species for desired chemi- 
cal reactions is a main quest in the field of heterogeneous 
catalysis. 1 Whether a species is active or not depends 
on many factors, including the nature of its interaction 
with the catalyst material, as well as the temperature and 
pressure at which the catalytic reactions take place. The 
main approach for gaining quantitative, microscopic un- 
derstanding of various properties and processes, has been 
to perform experiments under standard surface science 
conditions, i.e., under ultra high vacuum (UHV) and rel- 
atively low (room or lower) temperatures, with the hope 
that what is learnt may be relevant and/or extrapolated 
to the system under the "real" catalytic conditions of 
high temperatures and pressures. Theoretical studies of 
adsorption and reactions at surfaces have, at best, been 
carried out by first-principles approaches, which, while 
highly accurate, typically do not take into account the 
temperature or the gas-phase environment with which 
the surface is in contact. A recent focus, both experi- 
mentally and theoretically, is to attempt to study systems 
which more closely resemble those of true catalytic ones. 2 
Experimentally, this is done by performing in situ real- 
time measurements carried out under high temperature 



and pressure conditions, and by the systematic study of, 
e.g., well-characterized metal clusters deposited on sup- 
port mate ria ls a nd chemical reactions thereon (see e.g., 
Refs. Kl lll5ll()ll7lcSri . Such studies arc important in par- 
ticular for understanding those systems which exhibit a 
"pressure- and/or materials-gap" , i.e., where different be- 
havior is exhibited under low and high pressures, or when 
results for model single-metal catalysts do not represent 
those of the "real" catalyst, typically consisting of a dis- 
persed metal on a support material. Often the pressure- 
and materials-gaps are related and occur simultaneously. 

The oxygen-silver system is an example which ex- 
hibits a pronounced pressure-dependence: Despite its 
unique activity as the catalyst for ethylene epoxidation 
(C 2 H 4 + 5O2 ^ C2H4O, performed at 500-600 K and at- 
mospheric pressure) and partial oxidation of methanol to 
formaldehyde (CH 3 OH+±0 2 ^ CH 2 0+ H 2 - "oxyde- 
hydrogenation" , or by direct dehydrogenation, CH3OH 
^ CH2O+H2, performed at 700-900 K and atmospheric 
pressure), these reactions do not occur with any signif- 
icant probability under UHV conditionsi**i*'ii This be- 
havior of the oxygen-silver system might be related to 
a materials-gap, because for studies of the partial oxida- 
tion of methanol, drastic differences are observed in the 
performance for fresh silver catalysts and those which 
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have been exposed to high temperature (T > 873 K) and 
pressure (atmospheric) conditions. The temperature de- 
pendence of the conversion of the clean Ag material to 
the active catalyst exhibits a pronounced hysteresis, i.e., 
the conversion is different when raising the temperature 
from low to high, compared to from high to low. This 
was the case for high temperatures up to 873 K, as well 
as for lower temperatures up to 560 K and 750 K. These 
hysteresis effects were shown from scanning electron mi- 
croscopy (SEM) studies to be related to significant mor- 
phological changes which occur in the catalyst depending 
on the temperature, and that do not revert to the orig- 
inal structure when lowering the temperature, i.e., the 
changes are irreversible M> Structural changes in silver 
also occur when in contact with a pure O2 environment at 
high temperatures (^900 K) and atmospheric pressure, 
namely, facets with the Ag(lll) orientation resultii2ii£ 
This is presumably because under these conditions it is 
unlikely that there will be any oxygen on the surface and 
the (111) surface of Ag has the lowest surface energy. 

It has been suggested that the elusive sub-surface or 
"bulk-dissolved" oxygen species play an essential role in 
the above-mentioned catalytic reactions, which so far 
have escaped a precise characterizationiiSiiiiiiiSiiSiii In 
addition to this aspect, there are other controversial de- 
bates in the literature concerning the nature of the active 
species, e.g., atomic versus molecular (ozone- like) oxy- 
gen i 18 i 19 Another puzzling issue is how silver, as a noble 
metal, can function as a good catalyst when it only binds 
adparticles relatively weakly on the surface compared to 
the transition-metal catalysts to the left of it in the pe- 
riodic table. Catalysis clearly involves a number of com- 
plex processes, namely the dissociation of molecules and 
the creation of chemically active species, and subsequent 
interaction and reaction between the particles to form the 
product, which desorbs from the surface. Furthermore, 
it is possible (maybe likely) that some fragments of the 
reactant molecules (e.g., H or C) modify the catalyst ma- 
terial and thereby participate in the creation of the ac- 
tive species. In the present paper we tackle only part of 
the problem, namely, the study of the interaction of sil- 
ver with oxygen in order to identify and exclude possible 
active O species for the above-mentioned heterogeneous 
catalytic reactions. This is an essential prerequisite for 
modeling the full catalytic process. 

Although we study explicitly the (111) surface of sil- 
ver, our general findings are expected to be relevant for 
silver per se, and possibly to have implications for gold, 
also a noble metal oxidation catalyst^ which, for the 
(111) surface, exhibits a restructuring and chemisorption 
of oxygen atoms at elevated temperatures (500-800 K) 
and atmospheric pressure similar to Ag(lll). They 
may also relate to copper, which catalyses the oxidation 
of methanol to formaldehyde^ 

In our previous publications f2i2& the interaction be- 
tween oxygen and Ag(lll) was thoroughly studied by 
density- functional theory (DFT). In these works we 
considered almost all possibly relevant oxygen species, 



namely, pure on-surface chemisorbed oxygen, pure sub- 
surface and bulk-dissolved oxygen, surface- and bulk- 
substitutional adsorption, structures involving both oxy- 
gen on and under the top Ag layer, as well as surface- 
oxides and a molecular (ozone-like) species adsorbed at 
a surface vacancy. We studied these oxygen species for a 
wide range of coverages where the following results were 
obtained: At low coverages oxygen atoms adsorb on the 
surface. With increasing coverages, however, strong re- 
pulsive interactions develop between the partially nega- 
tively charged O atoms and the adatoms are driven to 
penetrate into the sub-surface region, and/or to give rise 
to a reconstructed, thin surface-oxide layer, when the 
coverage is larger than ^0.25 ML. For higher oxygen con- 
tents, thicker oxide-like films are predicted to form. 

We also found that the presence of sub-surface oxygen 
when bonded to the same surface Ag atoms as an oxygen 
atom on the surface, dramatically modifies the electronic 
properties of the on-surface oxygen and vice-versa. De- 
pending on the particular geometry, it can give rise to 
a notable stabilization or destabilization, but the ener- 
getically favorable structures induce a stabilization and 
a lowering of the density of states at the Fermi level. In 
contrast to the conclusions concerning the crucial role of 
bulk-dissolved oxygen based on interpretations of exper- 
imental results, we found that this species is energeti- 
cally unfavorable, as is bulk- and surface-substitutional 
oxygen, as well as the molecularly adsorbed ozone-like 
species at a vacancy, when the system is in thermal 
equilibrium^*^ On the other hand, the binding energy 
of atomic oxygen adsorbed at under-coordinated silver 
atoms is considerably stronger than on the terraces, but 
the concentration of such species obviously depends on 
the number of defects and the morphology of the sub- 
strate. 

Nevertheless, temperature and pressure can affect the 
stability of structures and this has not been investigated 
so far. Clearly this is important for gaining an under- 
standing of the system under catalytic conditions and 
this is the main topic of the present paper. Some of these 
results were reported briefly in Ref. |2|| Our study uses 
the approach of "a& initio, atomistic thermodynamics", 
that is, performing systematic DFT calculations for all 
O species that could conceivably be relevant, and take 
into account the effect of the environment, which we do 
via the pressure and temperature dep endence o f the oxy- 
gen chemical potential (see, e.g., Refs. l27128l29ll3(Mlll32r) . 
From the resulting free energies (obtained using our pre- 
viously obtained energetics—^), we derive the phase di- 
agram of oxygen at the silver (111) surface. Since our 
earlier studies2i2£ predict that oxide- like structures form 
for high O concentrations, we also investigate the sur- 
face and bulk properties of the most stable bulk oxide, 
namely, di-silver oxide, Ag20. 

The paper is organized as follows: The calculation pro- 
cedure, as well as the thermodynamical method are ex- 
plained in Sec. II, and in Sec. Ill results are presented 
for the bulk properties of Ag20. In Sec. IV, we investi- 
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gate various terminations of the Ag20(lll) surface as a 
function of the oxygen chemical potential. Section V de- 
scribes investigations into the transition from oxide-like 
structures (where the positions of the Ag atoms are com- 
mensurate with those of the Ag(lll) substrate) to the 
true oxide film which is laterally expanded in compari- 
son. In Sec. VI, the pressure-temperature phase diagram 
of the oxygen-silver system is presented and compared to 
experimental results. Finally, the conclusions are given 
in Sec. VII. 



II. CALCULATION METHOD 



phases. 

We consider a surface in contact with an oxygen atmo- 
sphere which is described by an oxygen pressure p and 
a temperature T. The environment then acts as a reser- 
voir, as it can give (or take) oxygen to (or from) the 
substrate without changing the temperature or pressure. 
We calculate the surface free energy, 



7(T,p) = (G - VA g MA g - Nono)/A 



(1) 



where No and iVXg are the number of oxygen and silver 
atoms, and A is surface area. The T and p dependence is 
mainly given by no, the oxygen chemical potential, i.e., 
by the O2 gas phase atmosphere: 



The DFT total-energy calculations are performed us- 
ing the pseudopotential plane wave method- 3 within 
the generalized gradient approximation (GGA)pi^ The 
pseudopotentials are generated by the scheme of Troul- 
lier and Martins with the same functionali^Siiil The wave 
functions are expanded in plane waves with an energy 
cutoff of 50 Ry. In the (lxl) surface unit cell, cor- 
responding to the periodicity of clean Ag(lll), 21 spe- 
cial k-points are used in the surface irreducible Brillouin 
zone (IBZ) for the Brillouin-zone integration; 38 Equiv- 
alent k-points to these are used for all of the surface 
structures studied for consistency, i.e., to maximize the 
accuracy when comparing the energetics of different cov- 
erages as calculated in different supercells. We employ a 
Fermi function with a temperature broadening parame- 
ter of ksT cX = 0.1 eV to improve the convergence, and 
the total energy is extrapolated to zero temperature, (fcs 
is the Boltzmann constant.) The structures are created 
on one side of the slab and the resulting induced dipole 
moment is taken into account by applying a dipole cor- 
rection.— For the adsorption structures on Ag(lll), we 
use five layers of Ag to model the surface and typically 
keep the bottom two or three layers fixed at the bulk- 
like positions and fully relax all the other atoms until 
the forces on the atoms are less than 0.015 eV/A . The 
oxide surfaces of Ag 2 0(lll) are simulated by symmet- 
rical slabs (i.e., with inversion symmetry). We tested 
thicknesses corresponding to five, seven, and nine metal 
layers, separated by 15 A of vacuum. Our results show 
that the systems can be safely described using five metal 
layers with inversion symmetry. 

In ab initio theory, the consideration of high tem- 
perature and high pressure can be achieved by explic- 
itly taking into account the surrounding gas phase in 
terms of ll ab initio, atomistic thermodynamics" (cf. e.g. 
Refs. I27ll2^l2g|l30ll3lll3^> . This is also an appropriate 
(first) approach to steady-state catalysis, which is of- 
ten run close to thermodynamic equilibrium (or a con- 
strained equilibrium) to prevent catalyst degradation. In 
the following we outline how the combination of thermo- 
dynamics and DFT can be applied to obtain the lowest- 
energy surface structures with a surrounding gas phase, 
thus enabling us to construct a (T, p)-diagram of the 
stability (or metastability) regions of different surface 



Ho(T,p) = l/2[E o °^ + f lO2 (T,p ) + fc B Tln^j],(2) 

where p° corresponds to atmospheric pressure and 
p,o 2 (T,p°) includes the contribution from rotations and 
vibrations of the molecule, as well as the ideal-gas en- 
tropy at 1 atmosphere. It can be calculated or taken 
from experimental values from thermodynamic tables - 
the difference between theory and experimental values is 
marginal in the temperature range of interest. The silver 
chemical potential is 



-bulk 



for adsorption on a silver substrate 



MAg I |(fl , A g2 0-buik - /^o),when bulk Ag 2 is present 

We use small symbols (s-Ag-bulk, 5A g2 o-buik) in Eq. 3 to 
indicate that these are free energies per Ag atom (top) 
and per Ag 2 unit (bottom). These two quantities are 
calculated by DFT-GGA. In Eq. 1, we have furthermore, 
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G , to obtain surface energies 



G 



LT 0/Ag(lll) 



G^tfiijN) to obtain adsorption (4) 



energies on Ag(lll) 



These quantities are also calculated by DFT-GGA. For 
the temperature dependence of p>o 2 (T,p°), we use ther- 
modynamic tables [Ref. |4Q| (i.e., experimental data), 
while the T = K value of Eq.^is 



Ho(T = 0K,p) = -E t 
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^Ag 2 0-bulk (5) 



In the above equation we use the experimental value of 
the heat of formation, #Ag 2 0-bulk' because of its very 
small value (0.323 eV at the standard state (room tem- 
perature and atmospheric pressure)^ or 0.325 eV ex- 
trapolated back to T — > OK), which implies that the error 
introduced in the later discussion by using the theoretical 
value (0.092 eV) would be significant. We note that the 
error in fljvgaO— bulk is due to the significant overbinding 
of 2 , and the underbinding obtained for bulk Ag and 
bulk Ag 2 0. 
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We choose the zero reference state of no(T,p) to be 
the total energy of oxygen in an isolated molecule, i.e., 
Mo(0K,p) = i^° 2 tal = 0. This will be called the 
"oxygen-rich" condition. The so-called "oxygen-poor" 
value of /io is taken as, 

poor j-itotal r) ptotal (r*\ 

MO ~~ £/ Ag 2 0-bulk — A ^Ag-bulk ■ \°) 

If the oxygen chemical potential becomes lower (more 
negative) than /^Q°° r , di-silver oxide starts to decompose. 
For oxides of different transition metals, the "oxygen 
poor" limit can vary significantly depending on the heat 
of formation. For example, what may be rather "oxygen 
rich" for RuC>2 (with Hf = 1.60 eV) is already too "oxy- 
gen poor" for Ag20 to be stable. The typical range of 
the allowed oxygen chemical potential is determined and 
limited by the pressures and temperatures used in indus- 
try and laboratories, namely, from several to thousand 
Kelvin, and from ultra high vacuum (< 10 -12 atm.) to 
several hundreds of atmospheres. 

For our investigations into various terminations of 
Ag20(lll) surfaces, we will also consider the energy re- 
quired to remove an oxygen atom from the surface and 
bring it to the vacuum. We call this the "removal en- 
ergy" , and it is defined as, 

^removal _ -E /Ag 2 0(Ul) + (E^ a O(lll) + ^o) , (?) 

where -Eo/Ag 2 o(iii) is the total energy of the surface ter- 
mination of interest, containing the O atom to be re- 
moved, -E , Ag 2 0(ni) * s t na * 01 the reference system, i.e., the 
surface structure remaining after the O atom has been re- 
moved. Eo is the total energy of a free O atom where 
the spin-polarization energy (—1.60 eV) is taken into ac- 
count. The removal energy reflects the bondstrength of 
the removed O atom to the surface. 

To help analyze the electronic structure, we also cal- 
culate the density of states, N(e), as well as the state- 
resolved DOS, or projected DOS, N a (e). The definitions 
are given in Ref. H2- 

III. BULK Ag 2 

The most stable silver oxide is Ag20, which decom- 
poses at 460 K under atmospheric pressured The crys- 
tallographic structure of di-silver oxide, Ag2 0, is shown 
m Fig.UJi. It has the cuprite structure with six atoms per 
cubic unit cell, (two oxygen atoms and four silver atoms) 
where the oxygen atoms form a body centered cubic lat- 
tice, while the metal atoms are located on the vertices of 
a tetrahedron around each oxygen atom, forming a face- 
centered-cubic latticed A characteristic feature, as seen 
in Fig. [Qi, is a linear oxygen-metal-oxygen coordination. 
The calculations for Ag 2 are performed with 35 spe- 
cial k-points in the irreducible part of the Brillouin zone 
(IBZ) of the six atom cell, corresponding to 1000 points 
in the full BZ. We used the same pseudopotentials and 
energy cut-off (50 Ry) as in our earlier works^ 2 ^ 



TABLE I: Calculated DFT-GGA and experimental bulk prop- 
erties of Ag2 0. The lattice constant, ao, the distance be- 
tween Ag atoms, dAg-Ag and oxygen atoms, do-o, and 
the bondlength between nearest neighbor oxygen and silver 
atoms, 6o-Ag, are given in A. The bulk modulus, B, (in 
Mbar), and the heat of formation, o-bulk' ( m e ^ P er 

oxygen atom) are listed as well. 

AO rfAg-Ag rfo-O fco-Ag B -gjgaO-bulk 

Exp. 4.72 a 3.34 a 4.09 a 2.04 a - 0.323 c 

4.74 b 3.35 b 4.10 b 2.05 b - 

Theory 4.91 3.47 4.25 2.13 0.647 0.092 

Reference b Reference H c Reference EH 



The calculated results, as well as available experimen- 
tal data, are presented in Tab. [I] The theoretical lat- 
tice constant of 4.91 A is 4% larger than the experi- 
mental one which is 4.74 A— or 4.72A42 The nearest- 
neighbor distance between Ag atoms in Ag2 is calcu- 
lated to be 3.47 A , which is larger than the (theoreti- 
cal) value in bulk fee silver which is 2.97 A . The calcu- 
lated O-Ag bondlength is slightly shorter than for oxygen 
chemisorbed on the Ag(lll) surface (compare 2.13 A to 
2.15-2.18 A , depending on the coverage), and similar to 
the O-Ag bondlengths of the (4 x 4) structure (depicted 
in Fig. |31 (bottom), discussed later and in Ref. l25|) which 
are 2.10-2.11 A for the upper O-Ag bonds and 2.14-2.16 
A for the lower O-Ag bonds (involving the surface Ag 
atoms). It is also similar to the O-Ag bond lengths of 
the favorable on-surface+sub-surface commensurate ox- 
ide structures (2.10-2.15 A) reported in Ref.|2^ [see e.g., 
the upper O-Ag-0 layer of Fig. |21 (top)] . 

The calculated heat of formation, obtained as the total- 
energy difference between di-silver oxide (Ag2 0) and the 
sum of the total energies of bulk silver (two Ag atoms) 
and the oxygen molecule (one O atom) is 0.092 eV. The 
difference to the experimental value of 0.231 eV is due to 
the significant overbinding of O2, which is 0.566 eV per 
oxygen atom, and the underbinding obtained for bulk Ag 
and bulk Ag2 O . This leads to some cancellation of errors 
but it is not complete. In table I, the calculated bulk 
modulus refers to the theoretical lattice constant. To the 
best of our knowledge, there are no theoretical values 
of the heat of formation with which to compare in the 
literature, nor lattice constants or bulk moduli. Both ex- 
periment and theory, nevertheless, show that formation 
of di-silver oxide is an exothermic process, and the small 
value is in-line with the low temperature of decomposi- 
tion noted above. The heat of formation of di-silver oxide 
is significantly smaller compared to the transition metal 
oxides for metals to the left of Ag in the periodic table, 
e.g. 1.60 eV, 1.19 eV, and 0.97 eV per oxygen atom (ex- 
perimental data at T — > K and atmospheric pressure) 
for Ru02, RI12O3 and PdO, respectively^ 

The total density of states is shown in the right panel 
of Fig. ^d. The peak located in the energy range —4.0 
to —2.0 eV with respect to Fermi energy arises mainly 
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FIG. 1: (a) The cuprite structure of bulk Ag20. The oxygen 
and silver atoms are indicated by the dark and pale grey cir- 
cles, respectively, (b) The band structure of Ag2 calculated 
at the theoretical lattice constant (left), and the correspond- 
ing density of states (DOS) is shown at the right. The energy 
zero corresponds to the Fermi level. 



panel of Fig.^J) at energy e = 0). We also performed cal- 
culations using the slightly smaller experimental lattice 
constant; this also did not yield a band gap. Using the 
augmented-spherical wave (ASW) method with an ex- 
tended basis set and the atomic sphere approximation 
(ASA) within the local-density approximation (LDA), 
Czyzyk et alm& calculated the projected density of states 
and band structure of Ag20, where a band gap 0.40 eV 
was found. However, no band gap was obtained by Deb 
and Chatterjee^i using the full-potential linearized aug- 
mented plane wave (FP-LAPW) method within the local- 
spin density approximation (LSD A). Their reported den- 
sity of states and band structure are very similar to those 
of our present work. In both of the two above-mentioned 
works, the lattice constant was not stated, but we believe 
the experimental one was used because our results cal- 
culated at the experimental lattice constant agree more 
closely. On increasing the energy cut-off and using harder 
pscudopotentials, namely, r s = = 2.20 Bohr and 
r p = 2.10 Bohr for silver and r s = r p = = 1.2 Bohr 
for oxygen, and an energy cut-off of 90 Ry, no significant 
change in the results was found (less than 0.02 eV for the 
heat of formation, 0.01 A for the lattice constant, and 
0.001 Mbar for the bulk modulus). Thus it appears that 
Ag2 is another semiconductor for which the LDA and 
GGA yield no band gap. Other materials for which this 
is the case are, for example, Ge^ InN^ and ScNiSi 

Calculations for Ag20 performed using the screened- 
exchange local-density approximation, 51 which affords a 
more accurate description of the band-gap compared to 
the LDA (or GGA), yield a band-gap of 0.73 eV which 
shows that this material is indeed a direct narrow band- 
gap semiconductor, in agreement with experiment. 



from the Ag-4d states, where there is only a small over- 
lap with the 0-2p band, which has main (but smaller) 
peaks in the regions of —6 to —4 eV and —2 to eV. 
The states in the former region are largely bonding and 
those in the latter are mainly antibonding in character, 
as determined by inspecting the spatial distribution of 
the wave- functions in these regions. The occupation of 
the antibonding states, and small overlap of O and Ag 
states, explains the low stability of the oxide. 

The band structure of Ag20 is shown in the left panel 
of Fig.^D. It can be seen that at the T point, the valence 
and conduction bands overlap, making Ag20 metallic, 
instead of semiconducting with a band gap of 1.3 eV 45 
as found in experiments. The spatial distribution of the 
wavefunctions in the conduction band close to the Fermi 
energy (investigated for several k-points), shows that the 
states are of antibonding nature, involving a hybridiza- 
tion of Ag-4d-0-2p- and Ag-4d-0-3s-like orbitals. The 
occupation around L of these bands could contribute to 
the underestimate of the heat of formation of Ag2 noted 
above. We note, however, that this down shift of the con- 
duction band has only a small effect on the DOS (cf. right 



IV. Ag 2 0(lll) TERMINATIONS 

In this section, we examine various surface termina- 
tions of Ag20(lll). This is interesting for two reasons, 
firstly in relation to learning more about oxide surfaces 
in general through comparison with existing studies, and 
secondly since our earlier works^ 4 *^ 5 . for oxygen adsorp- 
tion at the Ag(lll) surface predicted that for high oxygen 
concentrations, thick oxide-like structures form which 
contain half a monolayer of O between the Ag layers and 
quarter of a monolayer on the surface (see e.g., Fig. [3] 
upper right). On inspection of the geometry, we ob- 
serve that the structures correspond exactly to that of 
the Ag20(lll) surface, but where the positions of the 
Ag atoms are commensurate with those of the underly- 
ing Ag(lll) substrate, so the structure is laterally com- 
pressed compared to the true oxide [compare Figs. [21 and 
13 (upper)]. With respect to different terminations, addi- 
tional oxygen could adsorb on top of the Ag atoms. Such 
sites are called coordinatively unsaturated sites ("cus"), 
since compared to in the bulk, they are lacking their full 
O-coordination; for the present system, the cus Ag atoms 
are missing one O bond. Conversely, it could be favorable 
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FIG. 2: Top (upper) and perspective (lower) views of the 
atomic geometry of the Ag20(lll) surface. The silver atoms 
are represented by large, pale spheres and the oxygen atoms 
by small, dark spheres. The hollow and top sites of surface 
oxygen atoms are indicated, as is the "O-Ag-0 trilayer" which 
contains two O atoms and four Ag atoms per surface unit cell. 



that the surface has a metal-atom termination instead. 
We have therefore studied several different terminations, 
and calculated the surface free energy as a function of 
the oxygen chemical potential (cf. , Eq. [2} to determine 
which have the lowest energy. 

Figure [21 shows the two surface O atoms: The hollow 
site corresponds to the bulk stacking and the O atom in 
the top site makes this termination oxygen rich (without 
the O top atoms the site would be a "cus" site). The silver- 
terminated surface is that where both of these O atoms 
are missing. The stoichiometric surface can be thought 
of as a stacking of O-Ag-0 trilayers as indicated in Fig.|3 
which contains two O atoms and four Ag atoms per unit 
cell of the layer. We investigated the surfaces with oxy- 
gen in the hollow site, and the oxygen-terminated surface 
with oxygen atoms in the hollow and top sites. We also 
tested the silver terminated surface, as well as a termi- 
nation with oxygen atoms just in the top sites. 

The surface free energies of the considered surface ter- 
minations are shown in Fig. 0] The range of the chem- 
ical potential between the vertical lines corresponds to 
the (experimental) heat of formation of Ag20, which is 
0.325 eV, i.e. to the range in which bulk di-silver oxide 
is stable. It is clear that the surface with oxygen in the 
hollow site is the most stable one throughout the whole 
range of the chemical potential. The O terminated sur- 
face is very unfavorable. This is paralleled by the very 
large increase in the work function (by 2.84 eV with re- 
spect to the Ag-terminated surface, or 2.50 eV with re- 
spect to the stoichiometric surface with O in the hollow 




FIG. 3: Illustration of the atomic geometry of an oxide-like 
structure commensurate with the Ag(lll) surface. Top view 
(upper left) and perspective view (upper right) for an O con- 
centration of 1.25 ML. Ag atoms and O atoms are represented 
by large pale and small dark spheres, respectively. The struc- 
ture is very similar to the true oxide (Fig. but laterally 
compressed in comparison. Lower figure: sketch of the atomic 
structure proposed in the STM study of Ref.lH^for the (4 x 4) 
oxygen phase on Ag(lll). The surface unit cell is indicated. 
The oxygen atoms are represented by small dark circles, the 
uppermost Ag atoms by grey circles, and the intact plane of 
Ag(lll) atoms lying below the "O-Ag-O" trilayer are repre- 
sented as the white circles. Note that this surface structure 
is the same as a layer of Ag20(lll), except there is one Ag 
atom missing in the indicated cell, which would sit directly 
on top of the Ag atom in the (111) surface below (which is 
energetically unfavorable) . 



sites). Such a large work function increase indicates that 
the O atoms have a significant negative charge, and so 
there will be a considerable 0-0 repulsion. The surface 
terminated with only O atoms in the top site is even less 
favorable than the fully O-covered surface. 

The results described above are in contrast to those 
found for other (polar) oxide surfaces, e.g. RuO2(110) 32 
and a-Fe2O3(0001)j— where "full oxygen occupation" of 
available sites becomes favorable for high oxygen chem- 
ical potentials. The reason for this difference is because 
here the interaction between oxygen and silver is so weak 
that even extremely high values of the oxygen chemical 
potential cannot stabilize additional oxygen at the polar 
surface. 

In Tab. [H] various physical properties of the surfaces 
are listed. In particular, it can be seen that the removal 
energy of O-hollow is greatest when there is no additional 
oxygen at the neighboring top site (compare 3.97 eV with 
3.72 eV). Thus, the presence of O in the top site weakens 
the bond of the O-hollow atom to the surface, showing 
that there is a repulsive interaction between the O atoms 
in the hollow and top sites. We note that the value of 
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FIG. 4: Surface free energies, 7, as a function of the O 
chemical potential, for the various surface terminations of 
Ag20(lll): namely, Ag-terminated (dashed line), oxygen at 
the hollow site which corresponds to the stoichiometric sur- 
face (full line), and oxygen at both the hollow and top sites 
(dot-dashed line). The left and right vertical lines correspond 
to the "oxygen-poor- limit" and the "oxygen- rich- limit" , as 
explained in the text, and A/Lto is defined as 
using Eq.|^|for the latter. 
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3.97 eV is significantly stronger than 3.52 eV as obtained 
for O bonded to the clean Ag(lll) surface in the fee site 
at the same on-surface coverage of 0.25 ML as the O- 
hollow atoms shown in Fig. [21 It is, however, similar 
to the removal energies we found in our earlier work 25 
for the energetically favorable structures involving both 
on-surface and sub-surface oxygen and the compressed 
oxide-like structures, where the values are in the range 
3.90^4.01 eV. These bondstrengths can no longer be re- 
garded as weak, but are intermediate - and of similar 
order as those identified for the reactive species involved 
in certain chemical reactions^ 



A. Atomic and electronic structure 

The bondlength between oxygen in the hollow site and 
the surface silver atom for the energetically favorable 
Ag20(lll) structure is 2.11 A, which is slightly shorter 
than the corresponding bulk value (2.13 A) noted in 
Sec. III. It is similar to that of the energetically favorable 
compressed oxide-like structures (2.10-2.15 A ) [see Fig.0 
(top)] which have a very similar local geometry, but 
smaller than that for chemisorbed on-surface oxygen 
at Ag(lll) (2.15-2.18 A), as found in our previous 
worksi^iiS 5 . The bondlength between oxygen in the top 
site and the silver atom below of 1.96 A is significantly 
shorter because of its low coordination. Structurally, 
except for the Ag-terminated surface, where there is a 
large contraction of 10.3 % of the first metal interlayer 
spacing compared to that of the bulk truncated surface 
(calculated using the center of the layers), the relative 



TABLE II: Physical properties of the various surface termi- 
nations for Ag20(lll). £o™' ow is the removal energy of 
the oxygen atom which occupies the hollow site for the two 
systems considered (cf. Eq. 0. Adi2 and Atfe are the rel- 
ative changes in the first two metal interlayer spacings with 
respect to the bulk value (2.83 A), &o-A g is the bondlength 
of the surface oxygen atom to the uppermost metal atom, 
and $ is the work function. "o hollow " is the stoichiometric 
surface (O top removed, cf. Fig. and "O fu11 " denotes the 
structure shown in Fig. H where both O hollow and O top atoms 
are present. 



System 


£o"ow $ Ad 12 
(eV/atom) (eV) (%) 


Ad 23 
(%) 


60 -Ag 

(A) 


Ag-term. 


4.46 -10.3 


1.8 




Qhollow 


3.97 5.04 -0.6 


-0.7 


2.11 


QfaU 


3.72 7.30 -1.6 


0.4 


2.09 (hollow) 








1.96 (top) 



change in the metal interlayer spacings of the oxygen- 
terminated structures is modest, i.e., less than 2.0 %. 
For the stoichiometric surface, the rumpling of the 
silver atoms in surface metal layer is not negligible; the 
magnitude is of the order 0.10 A . This is due an outward 
movement of the Ag atom in the surface unit cell that is 
not bonded to the surface O atoms, i.e., the "cus" atom. 
These relaxations are smaller than those of the "true 
transition- metal" oxides (i.e., with incomplete filling of 
the d-band), e.g. a-Fe 2 O 3 (0001)2i and RuO 2 (110)^ 
which wc attribute to the weaker bonding between silver 
and oxygen. Similarly, for adsorption of oxygen on the 
Ag(lll) surface, we found only very small relaxations 
of the substrate metal layers^ in contrast to O on 
Ru(OOOl) 5 ^ where the O-metal bond is notably stronger 
and oxygen adsorption induces a significant expansion 
of the first metal interlayer spacing. The nature of 
the bonding of the di-silver oxide surface is similar to 
that of the bulk oxide: Oxygen and silver hybridize 
mainly via the Ag-4d- and 0-2p-like orbitals as shown 
in Fig. where the atom projected density of states of 
the first O-Ag-0 trilayer is shown. Compared to oxygen 
chemisorbed on the clean Ag(lll) surface (see Fig. 5 in 
Ref. 1^), the Ag-4(i band is significantly narrower, which 
is due to the increased distance between the silver atoms 
and the weak hybridization with 0-2p states. The PDOS 
in Fig. El is rather similar to that of the oxide-like struc- 
tures with same local bonding as presented in Fig. 12 of 
Ref. I24I In Fig. [S]two features, one at —5.53 eV and one 
at —0.66 eV, with respect to the Fermi energy, can be 
seen which are well separated from the center of the d- 
band. This is in contrast to the case of chemisorption on 
the metal substrate, where there is only one peak outside 
of the Ag-4d band, for example at —1.35 eV for oxy- 
gen in the on-surface fcc-hollow site at coverage 0.25 ML. 
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FIG. 5: Projected density of states (PDOS) of the first O- 
Ag-O trilayer (i.e., the stoichiometric surface): The surface 
oxygen (dashed line) which binds to three Ag atoms (it is also 
called hollow site oxygen). The Ag atom which is bonded to 
these O hollow atoms (solid line), and the sub-surface oxygen 
atom bonded to this Ag atom as well as to a single Ag atom 
in the layer beneath (dot-dashed line) [see Fig. |5] ignoring 
the O top atoms]. The energy zero is the Fermi energy, as 
indicated by vertical line. 



V. TRANSITION TO THE OXIDE 

Having studied the atomic and electronic structure of 
bulk Ag20 and determined the thermodynamically most 
favorable surface of Ag20(lll), we now wish to inves- 
tigate how a transition from the oxide-like structures, 
with positions of Ag atoms commensurate with those of 
Ag(lll), and laterally compressed by 14 % compared to 
the true oxide (as reported in Ref. |25|) , to the real oxide 
structure might occur. Such a commensurate oxide-like 
structure is depicted in Fig. [3] (top) for an O concentra- 
tion of 1.25 ML. In order to investigate such a transition, 
we study the difference of the total energy (per O atom) 
between oxide-like films (with commensurate Ag posi- 
tions to those of the metal substrate) and corresponding 
"real" (unstrained) oxide layers as a function of thick- 
ness. Because the structures can be thought of as the 
stacking of O-Ag-0 trilayers (see Fig. |2J), we calculate 
the energy difference with increasing numbers of trilay- 
ers. Then the number of layers at which the true oxide 
structure becomes energetically more favorable, defines 
a critical thickness for the transition to occur. In this 
approach, the coupling between the (strained) oxide-like 
film with the underlying metal substrate is not included. 
We note that this study may be somewhat "academic" 
since there is so far no evidence for such compressed oxide 
structures, however, this structure is a plausible inter- 
mediate and perhaps can even be stabilized under cer- 
tain experimental conditions. Furthermore, additional 
insight into the O/Ag system can be obtained by these 
calculations, and moreover, when compared to analogous 




Number of trilayers 

FIG. 6: Difference in total energy (per O atom) between a 
stack of O-Ag-O trilayers, commensurate with Ag(lll), and 
that of a corresponding oxide layer, as a function of the num- 
ber of trilayers in the stack (denoted by "+"). The commen- 
surate structures are just like the oxide, only laterally com- 
pressed by 14 % and vertically expanded by ~21 %. A nega- 
tive value indicates that the true oxide structure is energeti- 
cally more favorable. For comparison, analogous results (from 
Ref. l5a) are shown for Ru(0001) (denoted by "x") namely, the 
difference in total energy (per O atom) between a stack of 
O-Ru-0 trilayers [commensurate with Ru(0001)] and that of 
corresponding RuO2(110) layers as a function of the number 
of layers in the stack. The lines connecting the Ag and Ru 
results (full and dashed lines) are guides to the eye. Insets 
illustrate the atomic geometry of the bulk oxide structures 
where oxygen and metal atoms are represented by small dark 
and large pale circles, respectively. For Ru02, the O-Ru-0 
trilayer structure is also depicted (lower left), as well as its 
conversion, via an "accordian-like" expansion, into a "true" 
RuO 2 (H0) layer. 



studies performed for other O/transition- metal systems, 
in particular for O/Ru(0001)r^ various trends may be 
identified. 

The results are shown in Fig. A negative value in- 
dicates that the true oxide structures are energetically 
more favorable. We note that Fig. |H| assumes that the 
number of O atoms is fixed. However, in many (most) 
experiments that aim to study very high coverages, the 
actual thermodynamic variable may be no. 

For one O-Ag-0 trilayer, the "strained" [commensu- 
rate with Ag(lll)] oxide film is significantly more fa- 
vorable than the corresponding true oxide, as indicated 
by the relatively large positive value of 0.53 eV, but the 
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preference decreases quickly to 0.13 eV at three trilay- 
ers. On increasing the thickness further, the difference 
in the energy between the two structures decreases more 
slowly, and when the number of trilayers is nine, the true 
oxide film is energetically more favorable, as indicated 
by the negative value in Fig. Because each trilayer 
contributes half a monolayer of oxygen, the critical cov- 
erage for the transition from the commensurate oxide- like 
film to the true oxide film is around 3.5-4.5 ML. As the 
number of the trilayers approaches infinity, clearly the 
curve approaches the value corresponding to the energy 
difference (per O atom) between the bulk oxide and the 
"strained" [commensurate with Ag(lll)] bulk material. 

The electronic properties of the "strained" [commensu- 
rate with Ag(lll)] oxide film are very similar to those of 
the true oxide as can be seen from comparison of Fig. 12 
in Ref. with Fig. [5] This at first may be somewhat 
surprising since the lateral dimensions of the commen- 
surate, oxide-like structures are smaller, by 14 %. We 
find, however, that much of the associated strain energy 
is relieved by expansion of the Ag layers perpendicular 
to the surface to accommodate the O atoms in between. 
Through this relaxation, the lattice volume per atom is 
92 % of that of the bulk oxide. 

In Fig. H3 analogous results for the phase transition 
from O-Ru-0 trilayer structures (commensurate with 
Ru(0001)) to the corresponding real RuO2(110)^ oxide 
layers are also presented for comparison. Bulk ruthenium 
dioxide has the rutile structure. The atomic geometry 
of RuC>2(110) is shown as an inset in Fig. as is that 
of the (commensurate) trilayer and its expansion into a 
corresponding true RuC>2(110) layer (bottom, left). The 
commensurate trilayer structures are, laterally, 15 % and 
35 % smaller than for RuO2(110) (i.e., in the a and b 
directions of the rectangular (110) surface unit cell). In 
this case the transformation from commensurate trilayer 
to corresponding oxide layer, involves not only a lateral 
expansion like for the di-silver oxide structures, but also 
a "rotation" of the O atoms about the Ru atoms, as seen 
from the inset (explained in more detail in Ref. l55f) . In 
particular, every second O-Ru-0 component rotates to 
be horizontal to the surface, while the alternative compo- 
nents rotate so that they are perpendicular to the surface. 
In contrast to the Ag-0 system, the phase transition into 
RuC>2(110) oxide happens after only one O-Ru-0 trilayer. 
Also, for Ru(0001), each trilayer contributes two mono- 
layers of oxygen, so the critical coverage is around 2- 
4 ML. Because of the considerable energy gain due to 
formation of the real Ru02 oxide layers, compared to 
continuous of stacking of the O-Ru-0 trilayers, the for- 
mation of the oxide structure will certainly take place 
once it is kinetically possible. 

It is interesting to note that the energy differences of 
oxygen in the trilayers and the corresponding bulk metal 
oxide layers, exhibits a rather similar dependence on the 
number of layers for both O / Ag and O /Ru, except shifted 
in energy. This makes sense since we consider the en- 
ergy per O atom, and as the stack becomes larger, there 



are more and more O atoms in the center of the stack 
which have atomic environments which are similar (and 
thus similar energies), i.e., the affect on the energy due 
to the surfaces of the layers, diminishes. In view of this 
similarity, together with the fact that the heat of for- 
mation of e.g., the Rh (1.19 eV) and Pd (0.97 eV) ox- 
ides lie between the extremes of Ru (1.60 eV) and Ag 
(0.33 eV) oxides, we may expect that investigations of 
analogous transitions from commensurate structures to 
the bulk oxide phase for Rh and Pd would exhibit a be- 
havior intermediate to those of Ru and Ag. This would 
imply the trend that the transition from intermediate, 
commensurate, two-dimensional oxide films to the true 
oxide phase, occurs sooner (i.e., for thinner films) for the 
elements more to the left in the periodic table, like Ru, 
with larger heats of formation, compared to those to- 
wards the right, which have smaller heats of formation 
like Ag. 

From a recent trend study on the oxidation at the basal 
plane of late Ad transition metals (from middle to right in 
the periodic table, i.e., Ru, Rh, Pd, Ag)^& it was found 
that oxygen starts to enter the sub-surface region only 
after a certain critical coverage has been reached on the 
surface. Namely, at coverages beyond which the heat 
of oxygen adsorption on the surface becomes unfavorable 
compared to the heat of formation of the bulk oxide. The 
reported critical coverages are about 1.0 ML, 1.0 ML, 
and 0.75 ML, and 0.25 ML for Ru, Rh, Pd, and Ag, 
respectively. Beyond these critical coverages, onset of 
bulk oxide formation was proposed to take place. Clearly, 
however, due to kinetic hindering, bulk oxide formation 
may not be able to occur "immediately" and then it may 
proceed via the formation of oxide-like layers as described 
above. 



VI. THERMODYNAMICS OF THE 
OXYGEN-SILVER(lll) SYSTEM 

We now turn to investigate the affect of the temper- 
ature and pressure on the stability of the various oxy- 
gen/silver structures, the energetics and atomic geome- 
tries of which were reported in our previous publica- 
tions^^ in particular, the energies are summarized in 
Fig. 14 of Ref. |M 

To incorporate the effect of temperature and pressure, 
we calculate the Gibbs free energy as a function of the 
oxygen chemical potential (cf., Eqs. ^ and |2J similarly 
to the oxide surfaces, except here the reference system 
is the clean metal substrate. As discussed above, the 
"oxygen poor" limit is that below which di-silver oxide 
is unstable. For oxygen adsorption at the Ag(lll) sur- 
face, however, this limitation doesn't necessarily exist. 
The results are shown in Fig. [7] where we only show the 
low energy structures. The dependence of the energy on 
fio, namely the slope of the line, reflects the O cover- 
age. The higher the coverage, the steeper the gradient. 
We also include the corresponding temperatures for two 
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FIG. 7: Surface free energies for O at Ag(lll) for the various 
low energy structures as a function of the O chemical poten- 
tial, where Afio is defined as fio — l/2^° tal , using Eq. 
for the latter. On-surface chemisorbed oxygen with cover- 
ages 1/16 ML, 1/9 ML, and 1/4 ML (labeled by these cover- 
ages), the (4 x 4) structure (accordingly labeled "0.375 ML"), 
commensurate oxide-like structures with coverages 0.75 ML, 
1.25 ML, and 2.25 ML (as labeled by the coverages). The 
corresponding temperatures are given for two selected pres- 
sures (cf. Eq. 13, one corresponding to UHV conditions and 
the other to the atmospheric conditions of catalysis. 

given pressures, namely, that which corresponds to UHV 
and the other characteristic of the atmospheric pressure 
employed in catalysis. 

It can be seen that for low values of the oxygen chemi- 
cal potential, i.e. < —0.78 eV, the clean Ag(lll) surface 
is thermodynamically most stable. For higher values, 
(i.e., in the range of —0.78 eV to —0.63 eV), a low cover- 
age (1/16 ML) of oxygen on the surface in fcc-hollow sites 
is thermodynamically most stable, and for values greater 
than —0.63 eV and less than —0.54 eV, slightly higher 
coverages of adsorbed O (1/9 ML and 1/4 ML) become 
most stable. Increasing the oxygen chemical potential to 
values greater than —0.54 eV leads to the formation of the 
reconstructed (4 x 4) phase (labeled "0.375 ML"), which 
is stable up to A^o = —0.325 eV, i.e. to the "O-poor" 
limit (left vertical dashed line), where the bulk oxide be- 
comes stable. The atomic geometry of the (4 x 4) phase 
is illustrated in Fig. 01 (bottom). The atomic structure 
is very similar to a trilayer of Ag20(lll), except with 
one Ag in the (4 x 4) cell missing. The basis for this 
proposed structure (cf. Refs. 152115715 8 59 60) is the ob- 
servation that the area of the (\/3 x v / 3)i?30°-Ag 2 O(lll) 



surface unit cell is practically four times that of clean 
(1 x l)-Ag(lll); the difference in the length of the lattice 
vectors is only 1.2% (0.3% from the experimental lattice 
constants). Although the mentioned surface unit cells 
are (practically) commensurate, the positions of the Ag 
atoms within the x a/3)-R30° oxide surface cell are, 
all except one, incommensurate with those of the Ag(lll) 
surface. 

It can be seen that for values of A/io greater than 
—0.54 eV, the commensurate oxide-like structures (see 
e.g., Fig. |3J| (top) are also quite favorable (labeled by 
the O coverages, 0.75, 1.25, 2.25 ML). It can be noticed 
that the thicker the film, the more stable the structure 
is (i.e., the steeper the curve). Clearly, in the limit of 
an infinitely thick film, the corresponding line would be 
very close to the left vertical one. This is so because 
the average O binding energy in the thick commensu- 
rate oxide-like structures is very similar to that of O in 
the bulk oxide as described above. From Fig. it can 
be concluded that the bulk oxide is the thermodynami- 
cally most favorable structure once the oxygen chemical 
potential is higher than the "O-poor" condition, which 
is expected. However, the kinetics of oxide formation 
can be very slow and the thickness of the oxide film will 
be limited by O and metal bulk diffusion. Thus, it is 
unlikely that full thermal equilibrium of O2 with bulk 
di-silver oxide will be achieved. 



A. The O-Ag(lll) phase-diagram 

The results discussed above can be displayed as a con- 
tinuous function of temperature and pressure by con- 
structing a two-dimensional phase diagram, if we only 
consider the lowest energy structures. This has been 
done in Ref . |26| and the results are shown again in Fig. [S] 
Before giving a discussion, it is necessary to consider the 
magnitude of the error bars in the temperature and pres- 
sure. An indication can be obtained from knowledge of 
the stability of the bulk oxide, i.e., that it decomposes 
at 460 K under atmospheric pressure. From Figs. and 
OH the calculations predict that the bulk oxide decom- 
poses at ^350 K at atmospheric pressure, i.e., around 
110 K lower than experiment. At 460 K, the bulk oxide 
is unstable when the pressure is lower than 10 3 atm., i.e., 
compared to 1 atm. as observed experimentally. The rea- 
son for this disparity may be due to systematic errors of 
the DFT-GGA approach and/or to neglect of the entropy 
contribution. In particular, we note that the experimen- 
tal —TS contribution for bulk silver at 450 K and 1 atm. 
is —0.48 eV per two silver atoms, and —0.70 for Ag 2 
per formula unit4i Thus, the entropy contribution from 
di-silver oxide, compared to silver bulk (the entropy of 
O2 has already been included in the chemical potential), 
is greater by 0.22 eV per two silver atoms, which will 
shift the left vertical line in Fig. [7| further to the left. 
Therefore, in Figs. and El we generally expect that the 
temperature is underestimated, but the error less than 
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FIG. 8: Calculated phase-diagram for the oxygen- Ag(lll) 
system. On the horizontal axis is the the temperature in 
Kelvin, and on the vertical axis is the pressure. 



110 K, while correspondingly, the pressure is overesti- 
mated, by less than three orders of magnitude. For the 
temperature values quoted below, we give an estimate of 
the upper limit in brackets, that is, 110 K higher than 
that obtained from Fig. |SJ In spite of these uncertain- 
ties, we believe that the relative stability of the various 
systems and our general understanding and conclusions 
will not be affected. 

Typically, vibrational contributions to differences in 
the Gibbs free energies of extended systems exhibit 
some cancellation^ However, when there are additional 
atomic or molecular species which are not present in the 
reference system, the situation may be different as there 
will be no effective cancellation of the contributions of 
such species between surface and reference systems (as is 
the case using the clean Ag(lll) as reference). We inves- 
tigate this aspect in more detail for the present system; in 
particular, we consider low coverages of on-surface oxy- 
gen and thin oxide-like structures. We calculate the dif- 
ference in the vibrational and entropic contributions to 



the Gibbs free energy, F VIb = E vlb - TS V 



as, 
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where (3 = l/k B T (see e.g. Ref. I28l32j) . The vibrations 
which enter Eq. [5] are those of the adsorbed O atoms. 
The calculated perpendicular vibrational frequency for 
0.25 ML of oxygen in fee sites on the surface is 50 meV 
(400 cm -1 ), and we assume that the lateral vibrations are 
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FIG. 9: Vibrational contribution (oxygen only, see text) to 
the Gibbs free energy, _F vlb , for various adsorbate structures 
on Ag(lll). i) The structure with a coverage of 0.25 ML of 
oxygen in fee sites on the surface, ii) for the structure with 
one O atom in the fee on-surface site and one O atom in the 
sub-surface tetrahedral site (labeled "oxide-like"), and iii) an 
estimate for the (4 x 4) phase (dot-dashed line). 



the same (this is a somewhat rough approximation but 
suffices for the present discussion, see below). Thus, in 
calculating A_F vlb (T, ui) we multiply by a factor of three 
corresponding to these three modes. For the thin com- 
mensurate oxide-like structure containing one O atom on 
the surface in fee sites and one sub-surface O in tetrahe- 
dral sites (i.e., a total coverage of 0.5 ML), the calcu- 
lated perpendicular vibrational frequencies are 56 meV 
(455 cm -1 ) and 65 meV (526 cm -1 ), respectively, for 
these two species. Again, assuming each O atom has 
three modes equal to the perpendicular ones, we obtain 
the results presented in Fig. [5] 

It can be seen from Fig. that for on-surface O, for 
temperatures up to about 800 K, the contribution is a 
maximum of about ±2 meV/A 2 . For lower coverages of 
on-surface O such as the 1/16 and 1/9 ML structures, the 
contribution will be smaller (assuming the frequency is 
similar). For the oxide-like structure at coverage 0.5 ML, 
the effect is larger, namely 5-6 meV/A 2 for tempera- 
tures up to about 400 K, thereafter decreasing. For the 
(4 x 4) thin oxide structure with coverage 0.375 ML, we 
expect the vibrational frequencies of the O species will 
be similar to those of the structure at 0.5 ML since the 
local coordination of the O atoms is similar, as well as 
the electronic structure. In addition, the vibrational fre- 
quency of O in bulk Ag 2 is reported to be 66 meV 
(530 cm _1 )jSi which is very close to that of the O atom 
in the sub-surface tetrahedral site of the 0.5 ML cover- 
age structure (65 meV) . Because of the lower O coverage 
of the (4 x 4) structure, the contribution to the free en- 
ergy will be smaller than the 0.5 ML structure by about 
0.375/0.5 =0.75, giving a maximum contribution at low 
temperatures of 4.5 meV/A 2 (dot-dashed line in Fig.|5J). 

For the (4 x 4) structure, in the temperature range 



12 



400 K to 730 K, it gives rise to a destabilization of less 
than 3.5 meV/A 2 (the value at 400 K, which decreases 
with temperature), which on inspecting the energy scale 
of Fig. [7| can be seen to be non-negligible. The relative 
energy difference, however between the on-surface species 
and (4 x 4) structure, which determines the sequence of 
the stable phases, can be seen from Fig. [5] to be affected 
by less than 1.9 meV/A 2 (1.9 = 3.5 - 1.6 meV/A 2 ) for 
temperatures greater than 400 K, which will therefore 
have a small effect. We note that if the lateral vibrations 
would be 30 % greater than the perpendicular ones for 
the on-surface and oxide-like structures considered, this 
would result only in a small increase in the contribution 
to the free energy, i.e. of 0.5 and 1.2 meV/A 2 , respec- 
tively, at T close to zero. 

The above analysis may be regarded as qualitative, but 
it gives an estimate of the magnitude of the changes in 
the free energy due to consideration of vibrational and 
entropic effects of the extended systems for the energeti- 
cally most important O/Ag structures. Clearly, for other 
systems with adsorbed species that have high frequencies 
such as OH (of the order 400 meV)^ the contribution 
to the free energy from vibrational and entropic effects 
will be significantly greater and will have to be taken into 
account in determination of the phase diagram. 

Summarizing the results contained in Fig. [H] we ob- 
tain the following predictions: Under UHV conditions 
(p = 10~ 12 atm.), chemisorption of oxygen on Ag(lll) 
with low coverage (less than 0.0625 ML) is stable up to 
about 360 K (incl. corrections: 470 K). This is in-line 
with experimental results^ where it has been found that 
for temperatures less than 490 K, 0.03 ML of oxygen 
can be adsorbed on the Ag(lll) surface. Under UHV 
conditions, the calculations also predict the stability of 
the (4 x 4) phase at the lower temperatures of 160 K 
(incl. corrections: 270 K) to 260 K (incl. corrections: 
370 K). This is also consistent with scanning tunneling 
microscopy (STM) results obtained under UHV condi- 
tions at low (4 K) temperatures^ where atomic resolu- 
tion images of the (4 x 4) phase have been achieved. The 
(4 x 4) phase, however, cannot be directly prepared un- 
der these conditions as it is known from experiment that 
temperatures greater than 400 K are required in order 
for the necessary atomic rearrangement take place. It 
is therefore clear that it is the kinetics that prevent the 
reconstruction and formation of this phase at low pres- 
sures and temperatures. From Fig. [HI it can be seen that 
the (4 x 4) phase should develop into the oxide if the 
oxygen pressure is raised or if the temperature is lower 
than 160 K (incl. corrections: 270 K). Such transitions, 
however, may be prevented or hindered by kinetics (e.g. 
diffusion barriers) since extensive mass transport is re- 
quired to form the oxide. 

At atmospheric pressure (p = 1 atm.), typical of that 
used in industry, the bulk oxide is calculated to be only 
stable up to around 350 K (incl. corrections: 460 K). 
From this point of view, it can be concluded that for in- 
dustrial reactions which take place at considerably higher 



temperatures, di-silver oxide cannot be the active phase. 
The (4 x 4) surface-oxide phase, however, is stable in the 
temperature range of 350-530 K (incl. corrections: 460- 
640 K) which coincides with that at which epoxidation of 
ethylene occurs. We therefore propose that such atomic 
O species are actuating the reaction. With regard to the 
catalytic reactions that take place at even higher tem- 
peratures, it can be seen that only very low coverages of 
chemisorbed O are stable up to 720 K (incl. corrections: 
830 K); O atoms adsorbed at under-coordinated surface 
Ag atoms are up to 865 K (incl. corrections: 975 K) (see 
below) . Therefore these species could possibly play a role 
in these reactions. 

It is interesting to note that regardless of the pressure, 
Fig.[S]shows that the ordering of the stable phases is the 
same, so that the same phases should be observable at 
different pressures if the temperature is adjusted accord- 
ingly. This indicates a way to bridge the pressure gap, 
i.e., by appropriate choice of the pressure and tempera- 
tures to stay within the same phase. However, it is clear 
that sometimes this will not be possible within reasonable 
times because of kinetic hindering for dissociation and 
sticking of O2 (i.e., delivery of sufficient concentrations 
of atomic oxygen to the surface), as well as to kinetic 
limitations for adequate atomic rearrangement. 

Experimentally, it has been reported that a "strongly 
bound" oxygen species (named 7 ) forms after high tem- 
perature (~900 K) and atmospheric pressure exposure of 
polycrystalline silver, as well as Ag(lll), to an oxygen 
environment Temperature programmed desorption 
spectra show that 7 desorbs at around 900 K and it 
is proposed by the authors to be the active species for 
the dehydrogenation of methanoli 2 ^ where it diffuses 
to the surface from bulk sites. At the high temperature 
of ~900 K in an oxygen atmosphere, silver reconstructs 
to form many (111) facetsi2ii^ at the surface and it is 
speculated that this means that oxygen diffusion takes 
place by a so-called "substitutional diffusion mechanism" 
which the authors argue has higher barriers than usual 
interstitial diffusion, and to give rise to the O2 desorption 
peak at the high temperature of 900 K. From our calcula- 
tions, this picture is difficult to understand, in particular, 
how sufficient amounts of bulk-dissolved oxygen can be 
formed. From our theoretical study, we know that at 
lower coverage, oxygen will prefer to stay on the surface, 
instead of in the sub-surface or bulk region, and for high 
O concentrations, oxide-like structures are preferred^ 
Furthermore, our results showed that substitutional oxy- 
gen adsorption in the bulk (i.e., O on Ag sites) is ener- 
getically unfavorable. 25 



B. O adsorption at under-coordinated metal 
sites :Ag vacancies 

The importance of steps and vacancies as active sites 
has been well documented, from both experimental and 
theoretical points of view^Sa In our previous studies we 
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investigated the binding of oxygen atoms on surfaces with 
pre-existing vacancies^* We found the binding energy is 
stronger than on the (111) terraces, e.g., for 0.25 ML ad- 
sorbed in fee sites on a (2 x 2)/Ag(lll) vacancy array 
the binding energy is 3.83 eV while it is 3.52 eV on the 
perfect terrace, at the same coverage. A similar behav- 
ior can be expected for adsorption at under-coordinated 
sites at step edges, as has been shown theoretic ally for 
oxygen at other transition metals, see e.g. Refs. |53 and 
1631 To investigate the affect of pressure and temper- 
ature on the stability of such O species, we calculated 
the Gibbs free energy of formation. We find that under 
atmospheric pressure, the O atoms are stable to around 
770 K (incl. corrections: 880 K) on a (3 x 3)/Ag(lll) va- 
cancy array, and even to 865 K (incl. corrections: 975 K) 
on a (2 x 2)/Ag(lll) vacancy array. The existence and 
concentration of these species, however, clearly depends 
on the number of defect sites at the surface. In view of 
the above, it could be speculated that the TPD signal 
of "0 7 " is due not to substitutional diffusion of oxygen 
atoms from the bulk to the surface as proposed, but due 
to O atoms adsorbed at under-coordinated surface silver 
sites. 



C. Adsorption of O3 at an Ag vacancy 

As mentioned in the introduction, in the literature, 
the existence of surface Ag vacancies has been proposed 
to stabilize a molecular- like ozone (O3) species, which 
has been suggested to be active for the partial oxida- 
tion of ethylene i& As reported in Ref. |2{| we calculated 
the adsorption energy of such a species on (3 x 3) and 
(y/3 x y/S)R30° vacancy arrays, corresponding to oxygen 
coverages of 1/3 and 1 ML, respectively. We found that 
these species were energetically unfavorable compared to 
the commensurate oxide-like structures described above. 
If, however, we assume that Ag vacancies exist, i.e., ne- 
glecting the energy cost to create the vacancy, then the 
binding energy is the most favorable of all for coverage 
1 ML. We can estimate the concentration of such Ag va- 
cancies using our calculated values of the vacancy forma- 
tion energy, together with an Arrhenius type equation: 
For a (V3 x y/3)R30° vacancy array with a vacancy con- 
centration of 1/3 ML, the formation energy is 0.53 eV; it 
is 0.44 eV for a concentration of 0.25 ML, and 0.55 eV for 
a concentration of 0.11 ML. Considering a temperature 
of 550 K, corresponding to that around which epoxida- 
tion of ethylene takes place, this implies a concentration 
of thermally induced surface vacancies of less than about 
10~ 5 — 10" 4 ML, which is very small. Furthermore, cal- 
culation of the free energy shows that this species is less 
stable than the others. Our results thus do not support 
the molecular ozone-like species as playing an important 
role, but point to oxygen atoms of the thin (4x4) surface- 
oxide as being the main species actuating epoxidation of 
ethylene. 



D. Thick oxide formation at Ag(lll) 

As seen from Fig. our calculations predict that thick 
oxide-like films can form at Ag(lll), but only for low 
temperatures [e.g., <350 K (460 K incl. corrections) at 
atmospheric pressure] . The extremely low sticking coeffi- 
cient of O2 at Ag(lll)fi 4 i as well as the found low thermal 
stability of these structures, may hinder their formation. 
For example, at high temperatures at which diffusion of 
O (or Ag) atoms can effectively occur, the material may 
become unstable. Furthermore, when using NO2 to de- 
liver O atoms to the surface, the substrate has to be held 
at temperatures higher than the temperature at which 
di-silver oxide decomposes, in order to desorb the NO 
molecules. Using, however, low temperatures and atomic 
oxygen or ozone (which readily dissociates at the surface) 
instead of O2 (or NO2), we predict that thick oxide struc- 
tures could be observed experimentally - providing that 
diffusion barriers can be overcome. 

Interestingly, it was found that di-silver oxide (Ag20) 
can be prepared using a mixture of ozone and O2 or mi- 
crowave generated atomic oxygen, at room temperature 
on poly crystalline silver foil J^ ' 66 i 67 i 68 The oxygen gener- 
ated in these processes is quite active, for example, small 
amounts of ozone (5 mol.% O3 in O2) can generate the 
O partial pressure equivalent to that of a pure O2 at- 
mosphere at 10 +12 atm. at a temperature of 465 K. 66 
Actually, in these mixtures of reaction gas, Ag20 can 
be stable to temperatures in excess of 773 Ki££ In the 
studies of Waterhouse et. aZ.j££ it has been found that 
at room temperature, oxidation of polycrystalline silver 
foil is improved significantly after the sample has been 
"scratched" , which generates a high concentration of de- 
fects which act as the sites for the oxide film nucleation, 
highlighting the importance of defects. Also diffusion of 
oxygen and silver atoms is seemingly improved signifi- 
cantly due to the non-uniformity of the growth. 

We note that subsequent to completion of our studies, 
very recent experimental work reported thick oxide for- 
mation at Ag(lll) by using a flux of atomic oxygen^ 
thus supporting our findings. 

VII. CONCLUSION 

Through density-functional theory calculations, and 
incorporating the effect of the atmospheric environment, 
we obtained the pressure-temperature phase diagram for 
the oxygen/Ag(lll) system, which reveals important in- 
sights into this system and its function as an oxidation 
catalyst: The calculations predict that a thin surface- 
oxide structure is most stable in the temperature and 
pressure range of ethylene epoxidation and we propose 
it (and possibly other similar structures) contains the 
main O species actuating the catalysis. Low coverages of 
chemisorbed O become more stable for higher tempera- 
tures and could also play a role in other oxidation reac- 
tions. The only species stable for temperatures in excess 
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of 720 K (incl. corrections: 830 K), are O atoms adsorbed 
at under-coordinated surface Ag atoms. Thick bulk-like 
di-silver oxide can safely be ruled out as playing an im- 
portant role in the technical catalytic reactions due to the 
low thermal stability causing decomposition at temper- 
atures lower than those used in the reactors. They are, 
however, the most stable structures for strongly O-rich 
conditions and low temperatures (e.g., up to 330-460 K 
at atmospheric pressure). Interestingly, the electronic 
structure of thick oxide surfaces are very similar to that 
of the thin (4 x 4) surface oxide structure, which raises 
the question of whether di-silver oxide might be used as 
a low-temperature oxidation catalyst. 

A molecular ozone-like species adsorbed at a surface 
vacancy, as had been proposed in the literature as be- 
ing the active species for the epoxidation of ethylene, 
is energetically unfavorable, as is bulk-dissolved oxygen, 
which strongly questions the latters hitherto thought cru- 
cial role in, e.g., the partial oxidation of methanol to 
formaldehyde. 

We also studied the bulk and (111) surface properties 
of Ag 2 0, the most stable silver oxide. The interaction 
of oxygen and silver is via a weak hybridization between 
Ag-4d-5s and 0-2p orbitals, with substantial occupation 
of antibonding orbitals, which explains its low stability. 
Various terminations of the Ag 2 0(lll) surface were in- 
vestigated as a function of the oxygen chemical poten- 
tial and it is found that the surface with oxygen in the 
hollow site (i.e., the stoichiometric surface) is most fa- 
vorable for all of the range of the O chemical potential 



considered, while the Ag-terminated surface and surfaces 
with additional oxygen adsorbed on top of exposed "cus" 
surface Ag atoms are unfavorable. On investigating the 
energetics of a transition from commensurate oxide-like 
structures to the "real" (unstrained) oxide, we find that 
the transition becomes favorable only at nine O-Ag-0 
trilayers, which corresponds to ^4.5 ML oxygen. 

The approach of a ab initio atomistic thermodynam- 
ics", i.e., of calculating free energies for many possible 
atomic geometries and predicting the lowest energy 
structures in contact with species in a gas-phase environ- 
ment, is a promising approach to investigate a system 
under different pressure and temperature conditions, 
in particular, from those of typical theoretical surface 
science, right up to those of technical heterogeneous 
catalysis. The obtained phase diagram also shows that 
if one stays within the region of one phase or along a 
phase boundary, then a bridging of the pressure gap 
should be possible. Furthermore, it may be expected 
that the regions which lie close to the boundaries 
in the computed surface phase diagram may exhibit 
enhanced thermal fluctuations, i.e., here the dynamics 
of atomistic processes may be particularly strong, and 
this may be particularly relevant for catalysis under 
realistic conditions. Furthermore, we note however, that 
an understanding of a full catalytic cycle, requires a 
kinetic modeling that includes all reactant species and 
intermediates, which for complex systems such as those 
discussed in the present paper, is not yet possible. 
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